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We have measured the heat capacity of a mix- 
ture of He* and powdered cerium magnesium ni- 
trate (CMN) from 0.02°K to temperatures as high 
as 0.1°K at pressures of 0.12, 1.10, 3.4, 6.1, 

11.8, 21.1, and 28.0 atm and for pressures great- 
er than that of the minimum in the melting curve, 
Pmin, at external pressures of 30.7, 32.8, and 
34.8 atm. The specific heat seems to be linear 

in T below temperatures of about 0.03°K only at 
the lowest pressures, but at higher pressures it 
is possible that the linear temperature dependence 
of the specific heat characteristic of a Fermi liq- 
uid has not been developed by 0.02°K. Assuming 
a Fermi liquid behavior at 0.03°K, we obtain at 
zero pressure a ratio of the specific heat to tem- 
perature which is considerably greater than that 
deduced from earlier experiments? and predict - 
ed by Brueckner and Gammel,‘ but in agreement 
with the theory of Goldstein.® We also deduce a 
lower limit for the isobaric expansion coefficient. 
We find no discontinuity in specific heat down to 
0.02°K for the full pressure range in the liquid 
although the theories*®*’ of cooperative transitions 
in He® predict an increasing transition tempera- 
ture with effective mass and hence with pressure. 
We find that the liquid phase of He® at pressures 
above Pmin and at low temperatures has a small- 
&f specific heat and a reduced thermal boundary 
resistance relative to that at somewhat lower 
Pressures. We were able to measure the thermal 
boundary resistance at various pressures in a 
more restricted temperature range than the spe- 
tific heat and found that it decreased slowly as 
Pressure increased. Details of these measure- 
ments will be given in a subsequent publication. 


The measurements on the solid revealed a small 
specific heat increasing as the temperature de- 
creased below 0.036°K, indicating that the ex- 
change interaction is much smaller than that pre- 
dicted by Bernardes and Primakoff.® 

The cell containing the mixture of He* and CMN 
had the same geometry as that reported previous - 
ly,? except that there was 0.585 g of CMN in the 
cell and that the heater consisted of a bifilar lead 
of 0.002-inch diameter manganin wire 6 inches 
long, the ends of which were silver soldered to 
0.004-inch diameter copper leads. These joints 
were inside the cell. The Epibond 100A° plug, 
which closed the cell, was sealed by Epibond 121 
which is nonmagnetic.*° 

The magnetic susceptibility of the cell was care- 
fully measured in the 0.3-1°K temperature range 
both with and without CMN inside. We were able 
to calibrate the thermometer to within 2% and to 
deduce that magnetic effects caused an error of 
less than 1.5%. 

The pressure was applied to the He® using a 
stainless steel Toepler pump. The pressure was 
set using a dead-weight tester™ and monitored 
using a Helicoid gauge.’ The hydrostatic pres- 
sures due to mercury and oil in the pressure sys- 
tem were measured and corrections were applied 
to the dead-weight pressures. The pressure was 
held constant to within 1% during a run. 

The quantity of He® in the cell below 1°K was 
measured to be 0.0080+ 0.0002 mole at 0.12 atm. 
The He* content was less than 4 parts in 10° and 
the tritium decay rate in the cell was less than 
300 disintegrations/sec. 

The low temperatures were produced in essen- 
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FIG, 1. Ratio of specific heat, C, to gas constant R (=8.317 joules/mole K°) for He® under pressure. The right ep 
ordinate is displaced upward in order to display the data better. Above 29 atm, the pressures indicated were those with t 
in the external system. 3 giv’ 
latior 
; =2.1$ 
tially the same way as previously described.!»?° molar volume was assumed to be that of the liq- of C/ 
Temperatures down to 0.02°K were obtained when uid just before the filling tube plugged. This. the s 
the cell was thermally coupled directly to the re- assumption is probably not precise, but it does condi 
frigerator using copper “coil foils.”“* =Tempera- not lead to errors in excess of 5% in C/R. uid te 
tures down to 0.04°K were obtained when a lead The data of reference 1, plotted on Fig. 1, have# .1., 
switch was inserted between the cell and the re- been corrected for an incorrect temperature cal- Sums 
frigerator. ibration resulting from the magnetism”® of SC-13 speci 
The data for the heat capacity were taken by C/RI 
the method of reference 1. With the tight ther- Table I. Density ratios of He® near T=0°K. =2.71 
mal coupling, the best data were obtained around Bel 
0.03°K; and with the lead switch, the best data P - P s ese 
were obtained from 0.06 to 0.08°K. (atm) p/bo.i2 (atm) p/Po,t2 “a 
The ratio of the specific heat, Cc, to the gas con- 0.12 1.000 21.1 1.338 near 
stant, R, is displayed as a function of tempera- 
ture, T, for various pressures in Fig. 1. In re- 1.10 1.033 28.0 1,393 
ducing the heat capacity data to C/R, we assumed 3.4 1.094 30.7° 1.422 
the volume of the He® in the cell to be constant. 6.1 1.148 32.9? 1.437 
The number of moles in the cell was determined b wher 
from the molar volumes of Sherman and Edes- 43.8 5.56 8.8 wa’ ume, 
katy, * extended to 0K by extrapolating hp ~~ 0/0 o.49 is the ratio of the density of He’ at pressure 
sity data of Lee, Fairbank, and Walker*® to 0°K p to that at 0.12 atm near T=0°K, The molar volume 
and higher pressures. The molar volume ratios of He® at 0.12 atm is 36.73 cm®/mole. p/po,12 is used on 
used for the reduction of the data are given in to reduce the heat capacity data to specific heats. . 
Table I. For pressures greater than Pmjn, the bThis is the pressure at which the He® plug formed. = 
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flexible silver micropaint*® used to join the heater 
to the copper leads inside the cell. Since the 

game coil system and geometry were used in both 
reference 1 and this experiment, the temperature 
calibration in reference 1 was assumed to differ 
from the present one only by the ratio of the mas- 
ses of CMN. Both the corrected calibration con- 
stant of reference 1 and that in the present ex- 
periment agree within 5% with that calculated 
from the known geometry, galvanometer flux sen- 
sitivity, and properties of CMN. 

A lower limit on the values of C/RT and m*/m 
near T =0°K was obtained by drawing the best 
straight line between the data near 0.03°K and 
the origin. The results for various pressures 
are displayed on Fig. 2. There is no experimen- 
tal justification for this procedure at the higher 
pressures. At zero pressure we find that C/RT 
=2.78 (K°)“*+ 10% and m*/m =2.82+10%, using 


the formula 
me ma(ey) (1) 


where all symbols have their usual meaning.’ At 
these temperatures (Cp -Cy)/Cp is negligibly 
small. This value of m*/m is to be compared 
with the extrapolation from 0.085°K of reference 

3 giving m*/m=2.00+ 0.05, and with the extrapo- 
lation from 0.054°K of reference 2 giving m*/m 
=2.19+ 0.13. These are consistent with the curve 
of C/R vs T bending in the region of 0.03 -0.04°K, 
the same temperature above which the thermal 
conductivity began to deviate from the Fermi liq- 
uid temperature dependence.’” The theoretical 
value of m*/m of Brueckner and Gammel‘ is 1.84. 
Summing the spin and nonspin contributions to the 
specific heat calculated by Goldstein,® we find 
C/RT =2.74 (K°)™, which corresponds to m*/m 
=2.78. 

Below 25 atm, the limiting values of C/RT are 
linear in p and may be used to obtain the isobaric 
expansion coefficient. Using the formula, valid 
near T =0°K, 


oe, Rf 8 (I 
ade RT|, @) 


T 


where Ap = (1/V)(8v/8T)» and v is the molar vol- 
ume, we find at zero pressure 


ba, /8T =0.08 (K°)~?, (3) 


a8 a lower limit. Brueckner and Atkins” obtain 
0.076 (K°)~? for this quantity, using results which 
gave the effective mass incorrectly. Goldstein® 
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FIG. 2. Variation of lower limits of m*/m and C/RT 
with pressure. The points on this graph were obtained 
by taking the slope of the best straight line through the 
data near 0.03°K and the origin. Below 25 atm, the 
slopes of the curves are (d/dp)(m*/m) = 0.074 atm™ 
and (d/dp)(C /RT) = 0.036 (atm °K)~'. 


predicts 0.103 (K°)~*. Brewer and Daunt’® calcu- 
late from their data 0.105 (K°)~?+ 25%. Recently, 
Rives and Meyer” found Ap = -(0.12+0.02)T be- 
tween 0.04 and 0.08°K at 0.18 atm. 

We measured the specific heat and thermal 
boundary resistance of liquid He® at pressures 
above Pmin and at temperatures of 0.02°K and 
somewhat above by applying pressures of 30.7 
and 32.8 atm near 1°K and then cooling through 
the solid and back into the liquid. In this region 
both the specific heat and the thermal boundary 
resistance are anomalously low, and the specific 
heat had the opposite pressure dependence from 
that at somewhat lower pressures. During meas- 
urements taken when a plug formed in the filling 
tube at a pressure of 32.8 atm, the specific heat 
became very high at temperatures higher than 
0.0302 0.001°K, indicating the liquid-solid phase 
transition. For an external pressure of 34.8 atm, 
we found no evidence of two-phase equilibrium be- 
tween 0.017°K and 0.036°K, and hence, we con- 
clude that our specific heat measurements for 
these conditions were for the solid only. The re- 
sults of Sydoriak, Mills, and Grilly* indicate 
that, if the molar volume in the cell were the 
same as that in the liquid just before solidifica- 
tion commenced, it would be impossible to have 
only solid He* in the cell. Hence, we conclude 
that during the formation of the plug in the capil- 
lary leading to the cell, additional He*® enters the 
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cell. This effect can be observed in the external 
pressure system. 

The results of specific heat measurements of 
solid He® for an external pressure of 34.8 atm 
are displayed on Fig. 1. Although the data are 
not precise enough to determine the temperature 
dependence of the specific heat, the specific heat 
does decrease as T increases. The theory of Ber- 
nardes and Primakoff* predicts, for a body-cen- 
tered cubic crystal, that C/R =3J*/k?T?, where 
J/k is the exchange energy in temperature units. 
Assuming that the data fit a 1/T? law, we find an 
average value for |J/k| of 0.0015°K. This is 
about 10 times lower than the value of J/k 
(= -0.02°K) calculated in reference 8 for a face- 
centered cubic crystal. Since in the neighbor- 
hood of the melting curve and at low tempera- 
tures the entropy of the solid is unexpectedly 
high and the entropy of the liquid is unexpectedly 
low, it is not surprising that the predicted’ 
maximum in the melting pressure curve has not 
yet been observed. The rate of heat transfer be- 
tween solid He® and the walls is comparable to 
that of the liquid and the walls. 

The pressure of the minimum in the melting 
pressure curve was measured several times by 
observing the sudden increase in temperature of 
the CMN when the He’® pressure decreased slowly 
through Pmin, and the He* plug broke. We find 
Pmin=29.0+0.1 atm. This agrees, within exper- 
imental error, with Ppyjn = 28.914 0.02 atm ob- 
tained by Sydoriak, Mills, and Grilly™ and Pmin 
=29.1+0.1 atm obtained by Lee, Fairbank, and 
Walker.*® It is slightly lower than the Pmin 
=29.3+0.1 atm obtained by Baum, Brewer, Daunt, 
and Edwards.”* 

We wish to acknowledge helpful discussions 
with Professor John Bardeen. We are grateful 
for the advice given us by Dr. H. A. Reich and 
D. Fitchen in the construction of the pressure 
system. Dr. R. F. Nystrém very kindly meas- 
ured the tritium contamination in the He*. We 
wish to thank R. Sarwinski and W. Abel for help 
in the construction and running of the experiment, 
and J. Morr for assisting with the reduction of 
the data. 
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PROPAGATION OF SOUND IN He*T 


W. R. Abel, A. C. Anderson, and J. C. Wheatley* 
Department of Physics, University of Illinois, Urbana, Illinois 
(Received September 18, 1961) 


We have studied the propagation of sound in He® 
over a range of pressures extending from 1.4 psi 
to 513 psi. At low pressures we have obtained 
the sound attenuation coefficient and have inter- 
preted it in terms of a viscosity, the results be- 
ing in satisfactory agreement with experiments 
on other transport properties in He*®. In a pres- 
sure region somewhat above the minimum, Pypin, 
in the melting curve, we find a sharp break both 
in the velocity of propagation of sound and in the 
temperature dependence of the attenuation coef- 
ficient. 

The experimental cryogenic apparatus is simi 
lar to others which have been described else- 
where.’ The sound cell was sealed into a non- 
magnetic epoxy resin bulb with a powder of ceri- 
um magnesium nitrate which served as the ther- 
mometer. The sound cell consisted of an “X”- 
cut quartz crystal of nominal frequency 5 Mc/sec 
separated from an optically flat fused quartz re- 
flector by a fused quartz ring of thickness 0.1890 
ecm and bore 0.5 cm with parallel and optically 
flat faces. The ring had four radial grooves 
ground in one surface to allow entry of He*. Meas- 
urements of velocity were made using standard 
ultrasonic pulse techniques, no fewer than two 
echoes being used for this purpose, although in 
some cases as many as fifteen echoes were used. 
The attenuation coefficient was obtained by fitting 
the echo heights, h, to the equation h =h,e~%*, 
where a is the attenuation coefficient and x is the 
distance traversed by the sound. 

At a pressure of 3.38 psi, measurements of the 
attenuation coefficient were made at frequencies 
of both 5.13 Mc/sec and 14.92 Mc/sec. Near 1°K 
the ratio of the apparent attenuation coefficients 
was nearly the ratio of the two frequencies, char- 
acteristic of geometric attenuation. However, be- 
low 0.1°K the ratio of the attenuation coefficients 
was, within experimental error, given by the ra- 
tio of the squares of the frequencies. This result 
is characteristic of viscous attenuation, for which 
the attenuation coefficient is given by*® 


a =(817y? /3pc*)n, (1) 


where v is the frequency, p is the density, c is 
the sound velocity, and 7 is the viscosity. The 
experimental attenuation data are shown in Fig. 1. 
At 5.13 Mc/sec the geometrical attenuation coef - 
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FIG. 1. Ultrasonic attenuation coefficient as a func- 
tion of temperature in liquid He* at a pressure. of 3.38 
psi and at frequencies of 5.13 Mc/sec and 14.92 Mc/sec. 
The line through the experimental data at 5.13 Mc/sec 
for temperatures less than 0.06°K arbitrarily has been 
drawn with the T~? temperature dependence predicted by 
the Fermi liquid theory.° 


ficient is approximately 0.2 to 0.3 cm™. Below 
0.05°K to 0.06°K the attenuation coefficient obeys 
a T~™ law to within the accuracy of our measure- 
ments, the resulting viscosity being given by n 
=2.8x10~°T~ dyne sec (K°)?/cm?. Using the re- 
sult that 7 = (§v,?7,,) $0(m*/m), where v, is the 
Fermi velocity and m*/m is the ratio of quasi- 
particle mass to bare particle mass, to evaluate 
the relaxation time for viscosity, T,, from the 
Fermi liquid theory, we find tT) =17x10-™°T~* 
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Table I. Velocity of sound in liquid He’ at various 
pressures and, for pressures less than 426 psi, for T 
less than 0,1°K. 








Pressure Velocity Frequency 

(psi) (m/sec) (Mc/sec) 
1.42 184.0 5.13 
8. 06 192.1 5.13 
13.97 198.5 5.13 
16.25 201.0 5.13 
50.8 231.2 5.13 
89.9 260.2 5.13 
174.3 306.8 14.92 
314.6 361.8 14.92 
414.1 394.2 14.92 
423.9 396.8 14.92 
425.9 397.9 14,92 
433.6" 397.9 14.92 
483, 2% 408.6 14.92 
513. 07 416.9 14.92 





®The velocities at these pressures were obtained with 
T~1°K. 


sec (K°)?. This relaxation time is 2.7 times that 
found for the thermal conductivity.* Using the 
values of c, p,° and® m*/m to compute T,/rT, from 
the Fermi liquid theory, we find the theoretical 
estimate to be t)/Tp=1.5. 

At both 25 McSec and 15 Mc4ec, sound meas- 
urements were attempted under excellent condi- 
tions at 0.008°K, but no echoes were observed. 

The velocity of sound in liquid He*® at various 
pressures less than Ppmjn and at temperatures 
of 0.1°K and below is given in Table I. We esti- 
mate an error of +0.1% in the pressures and 
+0.3% in the velocities. The extrapolated veloci- 
ty at zero pressure is 182.4+0.6 m/sec, which is 
to be compared with 183.9+0.3 m/sec given by 
Laquer, Sydoriak, and Roberts.’ Also given in 
Table I are velocities in the liquid near 1°K for a 
few pressures greater than Pmin- These probab- 
ly may be compared meaningfully with those be- 
low 0.1°K, for we found that at 424 psi the veloci- 
ty did not change, within experimental error, be- 
tween 1°K and 0.1°K. 

At pressures greater than Pin, a plug of solid 
He® forms above the cell so that subsequent meas- 
urements are made at constant volume. In what 
follows, the pressures referred to are those in 
the liquid just before the plug is formed, and the 
frequency of the sound is 14.9 Mc/sec. The most 
striking features of the measurements in this re- 
gion are assembled in Table Il. For pressures 


300 


of 435 psi and below, we found velocities at low 
temperatures which were a few percent larger 
than those observed for pressures just less than 
Pmin- Moreover, the attenuation coefficient had 
both a magnitude and temperature dependence 
characteristic of the liquid at somewhat lower 
pressures. As the temperature increased, ech- 
oes were observed until the melting curve was 
reached, at which point we usually observed “dou- 
ble echoes” corresponding to propagation of sound 
between transmitter and reflector at both liquid 
and solid velocities, the latter being about 475 
+10 m/sec. Eventually, only the solid velocity 
was observed, followed by more “double echoes” 
and then the velocity in the high-temperature liq- 
uid. 

At pressures greater than 439 psi, the behavior 
of the liquid is remarkably different. At very low 
temperatures the sound velocity exceeds 500 m/ 
sec, and hence is greater than the velocity in the 
solid. This effect was demonstrated strikingly at 
a pressure of 483 psi (not listed in Table II) where 
we observed a “double-echo” pattern characteris- 
tic of two-phase equilibrium and corresponding to 
a velocity ratio of 14:13 = (510 m/sec:475 m/sec) 
in which the relative size of the high-velocity 
component steadily increased as the temperature 
decreased from 0.035°K to 0.016°K. In addition 
to having a sound velocity about 25% greater than 
that in the liquid at somewhat lower pressures, 
the He® in the high-pressure region is character- 
ized by an attenuation which increases, rather 


Table II. Characteristic sound parameters in He® 
at low temperatures and at pressures greater than the 
minimum in the melting curve. 








a b Attenuation 
Pressure Velocity Temperature coefficient 
(psi) (m/sec) (°K) (em~') 
431 409 0.030 0.9 
431 405 0.1 0.5 
431 402 Melting curve 
435 407 0.026 1.3 
435 404 0.1 0.5 
435 396 Melting curve “99 
439 503 0.030 1.0 
439 485 0.1 1.6 
454 510 0.022 aa 
454 491 0.1 2.8 





“These are the pressures at which the capillary 


plugged. 
Accuracy of velocity measurements is about 1%. 
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than decreases, with increasing temperature. 
Above the temperature region near 0.1°K we 
were unable to make measurements because of 
the high attenuation. Hence, as we increased the 
temperature above 0.1°K we saw practically no 
echoes until suddenly large echoes characteristic 
of the solid became visible, corresponding to the 
transition back into the solid. 

At a pressure of 513 psi we observed echoes in 
the solid only, the velocity being constant within 
experimental error from 1°K to 0.017°K. Accord- 
ing to the work of Sydoriak, Mills, and Grilly,® 
if the molar volume of the He® in the cell after 
the plug was formed remained constant at that for 
the liquid just before the plug formed, it would be 
impossible to have only solid He® in the cell for 
any of the pressures for which we made measure- 
ments. At an external pressure of 513 psi, and 
possibly at lower pressures also, we observed 
echo patterns characteristic of the solid only. 
Hence, we conclude that during the process of 
forming a phug in the tube leading to the cell, a 
substantial amount of He® passes into the cell. 
This conclusion is supported by work reported 
in reference 6, where under similar conditions 
the heat capacity characteristic of the solid only 
was observed. 

In the above measurements the pressure of the 
He® below the plug was not known and probably 
varied with temperature. Moreover, above 0.1°K 
the large attenuation made sound measurements 
impossible. However, the measurements do indi- 
cate that below 0.1°K and in a narrow pressure re- 
gion near 437 psi there is a drastic change in both 
the velocity and attenuation of sound. Moreover, 





it has been observed recently® that in this same 
region the specific heat is considerably lower and 
has the opposite pressure dependence from that in 
the liquid at pressures just less than Ppin- Hence, 
it appears likely that in this high-pressure region 
we are dealing with a different phase of He* from 
that which occurs at lower pressures. 

We are indebted to Professor John Bardeen for 
helpful discussions. We wish gratefully to ac- 
knowledge the assistance of E. N. Koch, M. Kuch- 
nir, W. Reese, G. L. Salinger, and R. J. Sarwin- 
ski, who aided us in assembling and running the 
experiments and in treating the experimental data. 
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SUPERCONDUCTIVITY OF NIOBIUM-MOLYBDENUM ALLOYS BELOW 1°K 







J. K. Hulm and R. D. Blaugher 
Westinghouse Research Laboratories, Pittsburgh, Pennsylvania 


and 


T. H. Geballe and B. T. Matthias 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received September 18, 1961) 


The alloys between group 5 and group 6 transi- 
tion elements are of great interest due to the 
availability of a continuous range of solid solu- 
tion between superconductor and nonsuperconduc- 
tor.’ In previous work the systems V-Cr, Nb-Mo, 
Ta-W, and Nb-W were studied down to 1°K and it 
was observed that the transition temperature of 
the superconducting component was lowered by the 
addition of the nonsuperconductor. The depres- 
sion of T, was found to be roughly proportional 
to concentration. 

It is of interest to know whether the observed 
linear relationship between 7, and composition 
extends below 1°K, particularly since such a re- 
lationship appears to be in conflict with the pre- 
diction of the Bardeen-Cooper-Schrieffer (BCS) 
theory.? The theory, with some simplifying as- 
sumptions,°® gives 

T = 0.855 6 exp[-1/N(0)V ], (1) 
where @ is the Debye temperature, N(0) is the 
density of states at the Fermi surface, and V is 
the (positive) interaction potential giving rise to 
the superconducting state. 

The niobium-molybdenum system was chosen 
for detailed study below 1°K. A number of compo- 
sitions in the range 35 to 45 atomic percent mo- 
lybdenum were prepared by levitation melting fol- 
lowed by high-temperature annealing (2500°C). 
The composition and impurity content were de- 
termined by chemical, gaseous, and spectro- 
scopic methods. Superconducting transition tem- 
peratures below 1°K were determined in a helium- 
3 cryostat using an ac method.* Sharp supercon- 
ducting transitions were obtained after improving 
the metallurgical techniques. Transition breadths 
ranging from 0.01°K up to 0.1°K were found for 
most samples. The resulting data are shown in 
Fig. 1 with the region below 1°K shown in detail 
in the insert. 

The dashed curve in Fig. 1 is calculated from 
Eq. (1) using N(0) and @ values interpolated linear - 
ly between the known values for pure niobium and 
pure molybdenum and assuming the value of V for 
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niobium remains constant.® While such a proce- 
dure is arbitrary, the long exponential approach 
to absolute zero is characteristic of the theory, 
not the procedure. The present experiments show 
an almost linear dependence of T, which extrapo- 
lates to zero at a finite composition, in conflict 
with the above application of the simplified theory, 
The possibility of this difference between theory 
and experiment has been emphasized by Pippard.° 
It is further illustrated in Fig. 2 where we replot 
the same data as Fig. 1 in the coordinate scheme 
suggested by Eq. (1). In this case the BCS plot is 
practically linear whereas the very steep descent 
of the experimental curve is well brought out. 
The linearly interpolated values of N(0) and é 
require that V in Eq. (1), or more precisely 
N(0)V/y where y is the electronic specific heat 
coefficient, be the parabolic function of concentra- 
tion shown in Fig. 1. A more rapidly varying N(0) 
than assumed would rotate the V curve counter- 
clockwise about the point for pure Nb and thus les- 
sen the rapid descent of V near 60% Nb. However, 
as the extrapolation in Fig. 1 indicates, 7, almost 
certainly goes to zero at an intermediate concen- 
tration,’ and therefore V must go to zero. Any 
deviation of the true N(0) from linearity should be 
due to the minimum in the density of states as Mo 
is approached and therefore would increase the 
curvature of the V curve as it approaches zero. 
It might be expected that a rapid change in V would 
be reflected in other ways; however, the known 
physical properties of the Nb- Mo system show no 
drastic changes in the pertinent concentration re- 
gion. Since V is the difference between the at- 
tractive phonon and repulsive Coulomb potentials, 
it possibly is a delicate function whose variation 
is not manifested in other ways. However, the 
mere fact that a large variation of V has to be as 
sumed is contrary to the fact that T, is a simple 
function of N(0) alone throughout the periodic sys- 
tem.® Such can be the case only if V is essentially 
constant, as indeed was assumed by Pines’ in ex- 
plaining the periodic occurrence of superconductiv- 
ity in the transition elements, or is a simple func- 
tion of N(0) itself. We are left with the dilemma 
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FIG. 1. Variation of the super- 
conductivity transition tempera- 
ture with composition in the Nb- 
Mo system. The dashed BCS 
curve shows the transition tem- 
perature calculated from Eq. (1) 
using a constant V. The N(0)V/y 
curve shows the variation in V 
which must be assumed to make 
Eq. (1) give the measured (and 
extrapolated) curve. 
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that such a V is not consistent with the present 
data. 

We would like to thank G. Comenetz for his 
levitation melting of the samples, G. W. Hull for 
help with the measurements, and P. Noziéres, 

J. C. Phillips, and P. Morel for their helpful com- 
ments. 
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PARAELASTICITY, A MECHANICAL ANALOG OF PARAMAGNETISM 
Werner K4anzig* 


General Electric Research Laboratory, Schenectady, New York 
(Received June 30, 1961; revised manuscript received September 1, 1961) 


It is known that certain radicals imbedded in a 
solid matrix of high symmetry can rotate even 
at very low temperatures.’ An interesting ques- 
tion to ask is: Is it possible to limit this rotation, 
i.e., to “align” the radicals, by a mechanical 
stress applied to the host crystal? For imper- 
fections that interact rather strongly with the 
matrix it has been shown by various techniques 
(internal friction, ? electron spin resonance’) that 
alignment can be achieved. However, the stress- 
alignable systems that have been reported so far 
freeze in at relatively high temperatures, and 
the kinetics of the hindered rotation is given by 
the conventional thermal activation theory (Ar- 
rhenius plot). 

In the present paper we report on a radical in 
alkali halides that has a different kinetic behavior. 
It still “rotates” at very low temperatures, and 
yet it can be aligned by modest mechanical stress- 
es. It is the O,” radical, which substitutes for a 
halide ion. A paramagnetic resonance study has 
shown that the internuclear axis of the O,~ mole- 
cule ions has six equilibrium orientations, name- 
ly, along the [110] axes of the crystal.* In the un- 
stressed crystal all six orientations are equally 
probable. If a mechanical stress is applied, the 


equilibrium orientations still correspond to [110] 
axes; however, the six [110] axes are no longer 
equally populated. 

This system lends itself very well to an accu- 
rate measurement of the orientational distribu- 
tion by means of paramagnetic resonance for the 
following reasons: (1) The anisotropy of the g 
factor is sufficiently large so that the paramag- 
netic resonance spectra corresponding to the six 
equilibrium orientations can be separated. (2) The 
O, molecule ion is very little affected by the 
stress field (except for its orientation). Therefore 
the intensity of the paramagnetic resonance aris- 
ing from the molecule ions oriented along a given 
[110] axis is a direct measure of their number. 
Thus we measured the stress alignment as a func- 
tion of compressional uniaxial stress along [111], 
[110], and [100], and as a function of temperature 
in the range from 1.5°K to 33°K for the host crys- 
tals KCl, KBr, and KI. 

For static stresses the distribution of the O,” 
radicals over the orientations, that become non- 
equivalent through the stress, was found to obey 
Boltzmann statistics. The energy differences oc- 
curring in the Boltzmann factors are proportional 
to the stress X. Figure 1 serves to illustrate these 
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FIG. 1. Alignment A=(N,-N,)/(N;+N,) of O,~ molecule ions in KBr for uniaxial stress X along [111], versus 
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X/T. The solid curve is the calculated curve tanh(aX/kT) for a =3.16<107-™ cm’, 
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Table I. Anisotropic part Ar ij, of the misfit, as de- 
termined from the eupertnentel Boltzmann factors. 








Units 1078 cm, 
Ar ’ ay ° Ar ’ 
xX yy Zz 
KCl -0.04 -0.40 +0.44 
KBr -0.15 -0. 26 +0.41 
KI -0.32 -0.19 +0.51 





facts for uniaxial stress along [111]. This stress 
gives rise to two energetically nonequivalent ori- 
entations, 1 and 2, so that the alignment A=(N,-N,)/ 
(N, +N.) is given by the well-known Brillouin func - 
tion for two energy levels. 

The interaction of the O,~ radical with the ex- 
ternal stress was analyzed in terms of the con- 
tinuum theory of lattice imperfections.* The O,~ 
radical was treated as an ellipsoidal inclusion 
in an elastic medium. The z axis of the ellipsoid 
was chosen along the internuclear axis, and the 
x axis along the lobes of the p function of the un- 
paired electron.* If one assumes that the elastic 
constants of the inclusion are close to those of the 
matrix, one obtains for the energy change of the 
inclusion due to an external stress® 

U=¢na*(S_ Ar +S Ar +S Ar). (1) 
xX XX yy YY 22 22 
Here a is the radius of the site (halide ion), Sj; 
is the component of the stress along the axis 7 
of the inclusion, and Ar}j; is the difference be- 
tween the radius of the free molecule ion along 
its axis j and the radius a of the halide ion. The 


Table Il. Constant B of the relation T=B/T. Units 
sec degree. 








B 
Stress axis ({111] (110) {100} 
KCl 21 +3 9 +3 11 +3 
KBr 4? 
KI 61 +5 76 +5 77 +2 





experimental results are internally consistent 
with Eq. (1) and permit the determination of the 
anisotropic part Ar,,’ of the misfit, Av;;’=Ar,. 

- §(Aryy+ Aryy+4rzz). The result is given in 
Table I. The numbers show that the molecule ion 
is cigar -shaped with the longest axis along z. 

The kinetics of the alignment process can be 
studied by monitoring the intensity of paramagnet- 
ic resonance lines upon sudden application or re- 
moval of the stress, since the paramagnetic re- 
laxation time 7, is at least one order of magnitude 
shorter than the reorientation times.” The time 
T* for alignment decreases with increasing 
stress, whereas the time 7 for recovery is inde- 
pendent of the magnitude of the stress, as is ex- 
emplified in Fig. 2. Below 4.2°K the time 7, as 
well as T* for moderate stresses, is inversely 
proportional to the temperature 7, i.e., T=B/T. 
Thus the system can be considered as a mechan- 
ical analog of paramagnetism, not only with re- 
gard to its static behavior (Fig. 1), but also with 
regard to its relaxation behavior at very low tem- 
peratures. The experimental values of B are giv- 
en in Table II. 
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FIG. 3. Example of tem- 
perature dependence of the 
product TT. KI stressed along 
[100]. 
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The anisotropy of the reorientation kinetics is 
(within experimental error) consistent with the 
assumption that the reorientation involves two 
basic processes, namely, a 90° rotation in a (100) 
plane with the probability p,, and a 60° rotation 
out of a (100) plane with the probability p,. The 
analysis of the rate equations, using the B values 
in Table II, yields: 


for KCl, D, << po, 
p,/T=(15+3) x 107° sec™ degree™, 
for KI, p,/T=(3.8+0.4) x 10~° sec™ degree", 


p,/ T=(2.2+ 0.2) x 107° sec degree™. 


The B/T law suggests that the reorientation of 
the O,~ radical involves emission or absorption of 
one phonon. If this is true, the ultrasonic absorp- 
tion versus stress should exhibit resonance peaks 
(“paraelastic resonance”). Inverted populations 
can be generated by reversing the stress so that 
mechanically pumped phonon masers might be- 
come feasible. 

Above 4.2°K the time t becomes much shorter 
than B/T (Fig. 3). Usual thermal activation seems 
to take over. 


A detailed account of the present investigations 
will be published elsewhere. The author grate- 
fully acknowledges the cooperation of G. D. Wat- 
kins, W. G. Johnston, T. G. Castner, and G. E. 
Nichols. Special thanks are due to M. H. Cohen 
for active help at important stages of the inter- 
pretation. 
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DOUBLE RESONANCE AND NUCLEAR COOLING IN AN ANTIFERROMAGNET* 


A. J. Heeger,t A. M. Portis, Dale T. Teaney,!! and Gerald wittT 
Department of Physics, University of California, Berkeley, California 
(Received September 18, 1961) 


We wish to report the observation of the Mn™® 
nuclear resonance by a double-resonance tech- 
nique in the antiferromagnet KMnF,. In addition, 
by observing the nuclear spin temperature as a 
function of power absorbed by the electronic anti- 
ferromagnetic resonance, we have discovered an 
apparent cooling of the nuclei via their interac- 
tions with the electron system. 

Examination of the antiferromagnetic resonance 
of KMnF, in the liquid helium temperature range 
shows that even at these low temperatures there 
is a strong temperature dependence of the reso- 
nance field as is shown for one of the absorption 
lines by the experimental points of Fig. 1. We 
have interpreted this shift to be the result of a 
strongly temperature-dependent anisotropy field 
arising from the hyperfine interaction of the anti- 
ferromagnetic system with the Mn® nuclear mo- 
ments. This hyperfine anisotropy field has the 
form 


* « (A/gu,dL). 
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FIG. 1. Theoretical fit of the temperature dependence 


of the resonance field. The mechanism is the tempera- 
ture dependence of the hyperfine anisotropy field. 





The hyperfine constant A/gu p has been evaluated 
experimentally by Ogawa! from the paramagnetic 
resonance of Mn** in isomorphous compounds. 
Substituting Ogawa’s value of 92 oe, we obtain 
= 8.0/T oe. 

H mn /T oe 
The expression for the resonance field of the elec- 
tronic antiferromagnetic resonance is 


+ 2H 2+ (w/ )?+5H H - 2H H e 
. [3 A, ti E ly 


H,=-§H EM 4, 


0 A 
The anisotropy field H, is responsible for the 
canting in this material* and is equal to 1300 oe 
as determined from measurements of the mag- 
nitude of the weak moment. The field H A, which 
is associated with the cubic and axial anisotropy,° 
and the exchange field, Hr, were determined as 
3.6 oe and 1.6x10° oe, respectively, by fitting 
theory and experiment. From the above expres- 
sion it may be seen that the shift in resonance 
field resulting from the nuclear magnetization is 
of the order of (Hg/Ho)HA,. Thus a one-oersted 
nuclear field may produce a shift in resonance 
field by nearly one kilo-oersted. The completely 
satisfactory agreement between the experimental 
observations and this theory is shown in Fig. 1. 

The strong dependence of the antiferromagnetic 
resonance field on nuclear spin temperature sug- 
gested the possibility of observing the Mn** nuclear 
resonance in antiferromagnetic KMnF, by a double- 
resonance technique. If enough rf power is sup- 
plied at the nuclear resonance frequency to sat- 
urate the nuclear magnetization, a shift of the 
antiferromagnetic resonance to its high-temper- 
ature position is expected. From the measure- 
ments of Ogawa the expected nuclear resonance 
frequency is 


Vo= (A/h){S) = 630 Mc/sec. 


The sample was a small cube of KMnF,, approxi- 
mately 1 mm on a side, at the center of a small 
rf coil, slightly over 1 mm in diameter. The sam- 
ple and coil were placed at the bottom of a rec- 
tangular microwave cavity resonating in a TE,,, 
mode. A shift in the position of the antiferromag- 
netic resonance with rf power could be observed 
over a 42-Mc/sec range in frequency centered 
around 628 Mc/sec, in good agreement with the 
expected central frequency. In Fig. 2 is shown a 
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quency. It is found that when the rf field exceeds Mn® nuclei themselves very strongly inhibit the the | 
a critical value the electronic resonance line shifts sublattice excitation in the vicinity of the nuclear acti 
abruptly to a new position characteristic of a high resonance as long as there is a substantial nuclear The 
nuclear temperature. On reducing the rf field we moment. This theory indicates that it should be poss 
find that the resonance line remains at its high- harder to produce an initial saturation of the nu- orde 
field position down to rf power levels well below clei against their dynamical resistance than to chi | 
those required for the initial shift. maintain the saturation when there is very little thou 
By pulsing the 628-Mc/sec rf field on and then sitin 
off, it was possible to observe the recovery of the are 
position of the antiferromagnetic resonance line — ' , ; : P aun 
after removing the rf power. Such a recovery aa 7 tion 
curve after nuclear saturation at 2.1°K is shown 800 7 an e 
in Fig. 3. The relaxation appears to be nearly 700+ J proc 
exponential with a characteristic relaxation time be c 
of 50 msec. 4 600 7 fror 
The very large shifts indicate that the nuclear $ The 
magnetization is nearly reduced to zero at the rf 8 500 : neti 
levels employed. We estimate the interaction of = has 
the Mn®* nuclei via spin waves‘ to yield a charac- ~~ 400} 4 the | 
teristic time T, of the order of 0.15 ysec. With = urat 
the observed 7, of about 50 msec, the rf field for Ww 
partial saturation should be approximately 2 oe. 300 + 4 of P 
This is actually about the level at which the nu- his 
clear saturation jumps to nearly its full value. As autl 
is shown in Fig. 2 the nuclear saturation may actu- 
ally be maintained at rf fields well below this val- 200 I j \ | | | 
ue. We have analyzed the excitation of the Mn™ 0 10 20 30 40 50 60 = 70 
nuclei in some detail in an effort to understand the Time, msec 
complex behavior observed. We find that because FIG. 3. Recovery curve showing the shift of the anti- 


of the large parallel susceptibility of a canted anti- ferromagnetic resonance at a given time after turning 
ferromagnet,® the Mn® nuclear resonance is ex- off the 628-Mc/sec rf power. 
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nuclear moment, in agreement with the observed 
hysteresis shown in Fig. 2. The frequency range 
over which such behavior may be expected depends 
on the nuclear magnetization and bears very little 
relation to the actual linewidth 1/7,. 

Experiments have also been performed® in which 
the antiferromagnetic resonance is driven at high 
microwave levels without any direct nuclear ex- 
citation. Under these conditions we also find a 
shift in the position of the antiferromagnetic reso- 
nance. This shift is clearly associated with a 
change in nuclear magnetization. The apparent 
nuclear temperature as a function of microwave 
power for several lattice temperatures is shown 
in Fig. 4. At 4.2°K the nuclear temperature re- 
mains constant up to the power at which spin- 
wave breakdown occurs.® At high power levels 
the increased nuclear temperature presumably 
reflects the general increase in the level of ex- 
citation of the magnetic system. At 2.1 and 1.8°K 
there is a quite unexpected reduction in the nuclear 
temperature before the onset of spin-wave heating. 
At these lower temperatures and low power levels 
the nuclei are evidently cooled through their inter- 
actions with the off-equilibrium spin-wave system. 
The cooling of nuclear spins may be related to the 
possible development of an Overhauser effect in 
ordered magnetic materials as discussed by Ogu- 
chi and Keffer.” It should be emphasized that al- 
though the present experiments are unusually sen- 
sitive to a change in nuclear temperature, they 
are not especially sensitive to the development of 
auniform nuclear polarization. Such a polariza- 
tion would shift the resonance in the same way as 
an external field. Full nuclear polarization would 
produce a shift of only a few oersteds and would 
be completely dominated by the shift resulting 
from even a slight change in nuclear temperature. 
The recovery of the position of the antiferromag- 
netic resonance following microwave saturation 
has been studied. The relaxation appears to be 
the same as that observed following direct rf sat- 
uration of the nuclei. 

We wish to acknowledge the helpful comments 
of Professor T. Nagamiya and the hospitality of 
his laboratory extended to one of us (A.M.P.). The 
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FIG. 4. Spin temperature as a function of microwave 
power absorbed by the sample. The spin temperature 
is determined by the resonance field. 
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OBSERVATION OF THE MAGNETICALLY INDUCED MAGNETOELECTRIC EFFECT 
AND EVIDENCE FOR ANTIFERROMAGNETIC DOMAINS 


G. T. Rado and V. J. Folen 
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In this Letter we report the first observation of 
the magnetically induced magnetoelectric effect 
[henceforth denoted by (ME) 7] in an antiferromag- 
netic material. When correlated with the previ- 
ously observed’? electrically induced magneto- 
electric effect [henceforth denoted by (ME); ], this 
observation provides a necessary verification of 
the existence of EH -type terms*** in the thermo- 
dynamic potential of a spin-ordered substance. 

We also report the discovery of magnetic anneal- 
ing effects in an antiferromagnet. In conjunction 
with other structure-sensitive phenomena which 
we found by means of magnetoelectric (ME) meas- 
urements, these magnetic annealing effects will 
be shown to provide the first strong evidence for 
the existence of antiferromagnetic “domains.” 
Previous work® on such “domains” was largely 
speculative. 

The thermodynamically derived* constitutive re- 
lations appropriate for the (ME)y effect are 


Di =€)E +a)H, D. = £ 


+aH 
lz lz x,y t+ x,y 


1 x,y’ 
(1) 


where a, /47 and a,/4n are magnetoelectric sus- 
ceptibilities and the z axis is chosen to be along 
the principal axis (c axis) of a uniaxial crystal. 

In our experiments on Cr,O, the applied electric 
field was zero. Using x-ray-oriented single-crys- 
tal disks whose cylindrical axes z’ were at an an- 
gle 6 (@=0, 90°, and 6, =30°) to the z axis, we 
switched on a magnetic field H along z’ and meas- 
ured the voltage V which was induced between the 
silvered plane surfaces of the disks. For fields 

up to the highest values used (H = 2500 oe), V was 
found to be proportional to H. Although V decays 
with a relaxation time t determined by the conduc - 
tivity and dielectric constant of Cr,O,, the value 
of t was found to be sufficiently large to permit 

us to measure V with a vacuum tube electrometer 
(Keithley Instruments, Inc., Model 610R). In the 
vicinity of room temperature the low-voltage meas- 
urements could not be made because of drift in 

the apparatus. The actual voltage V may not be 
equal to the open-circuit voltage V° (because the 
resistance of the disks may not be small in com- 
parison to the input impedance of the electrom- 
eter), but it is clear the V must be proportional 
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to V° and hence [because of Eq. (1)] proportional 
to a. 

Our experimental results show that the tempera- 
ture dependence of V, and V, is essentially the 
same as that of the a, and a, which we deduced 
from our (ME); effect? experiments and described 
satisfactorily in a recent theoretical® study. As 
shown in Fig. 1, the temperature dependence of 
1Vg. | and lag 4PP| [which were measured by 
means of the (ME)y and (ME), effects, respec- 
tively] are also in agreement. Since the values 
of antiferromagnetic permeabilities are close to 
unity, we find that ag@PP (the measured value of 
ag) is given by the relation ag®PP =[1 - (Nz /An)lag, 
where N,, is the axial demagnetizing factor of an 
oblate spheroid having the same axial ratio as the 
disks used in the ME measurements. Further- 
more, we observed that the signs of Vo, and 
ag app reverse at the same temperature. We 
also found that the three ratios IV t/a y PPI, 
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FIG. 1. Comparison of the temperature dependence of 
lag, @PP| and the induced voltage (|Vg,|). These quanti- 
ties were measured by means of the (ME). and (ME)y 
effects, respectively. Both ag, @PP and Vo, change sign 
near 154°K. The angle 4; is about 30°. 
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Iv,1/la,2PPI, and 1Vg,1/lag,4PPi are equal to 
each other provided the V and a@PP contained in 
a given ratio are measured on the same sample. 

Using both the (ME)y and the (ME)g effects, we 
observed that the signs and magnitudes of a and 
a,, but not their temperature dependence, may 
vary (a) from sample to sample, and (b) upon cool- 
ing a sample through the Néel temperature either 
with or without the presence of a magnetic field. 
In view of (a), it is understandable that the rela- 
tion ag =a 4 cos?@ +a, sin?@ (which follows from 
the form of the a tensor of Cr,O,) is not obeyed 
if the a’s are measured on different samples. 

The magnetic annealing mentioned in (b) caused 
the a’s to increase by factors as large as 3000. 
In addition, an almost complete “erasure” of a i 
was accomplished by cooling a z l|z’ sample 
through the Néel temperature in the presence of 
a 60-cycle/sec magnetic field. 

All these observations may be interpreted by 
postulating the existence of antiferromagnetic 
“domains.” As to the specific nature of the do- 
main structure, we tentatively suggest that in the 
case of Cr,O,, there are two kinds of domains 
and that they are characterized, respectively, 
by each spin in a domain pointing toward, or away 
from, the nearest oxygen plane which is perpen- 
dicular to the c axis. The difference between 
these two kinds of domains results from the asym- 
metry of the Cr*+* sites with respect to the oxy- 
gen planes. This asymmetry produces linear 
terms in the crystalline electric potential, and 
in a possible mechanism*® of the ME effects these 
terms play an essential role and determine the 
signs of the a’s. Consequently, the suggested 
domain structure explains the fact that the ob- 


served @’s can have either sign and that (due to 
partial cancellations resulting from the presence 
of two types of domains) they can have variations 
in their magnitudes. Since the two kinds of do- 
mains differ solely by a 180° reversal of all spins, 
it is clear that the intrinsic a@’s associated with 
each domain should have the same temperature 
dependence. Thus we have a natural explanation 
of the striking fact that the temperature depend- 
ence of each observed ag does not vary even 
though its sign and magnitude may vary. Finally, 
we note that the large increases in the magnitudes 
of the a’s which can be produced by magnetic an- 
nealing from a temperature just above the Néel 
point are difficult to explain without postulating 
antiferromagnetic domains. 

We wish to thank Miss J. W. Radue and Mr. 
M. J. Marrone for help in taking some of the 
data. 
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EXCITATION OF SPIN-WAVE RESONANCE BY MICROWAVE PHONONS 
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The interaction of acoustic vibrations with the 
spin waves in ferromagnetic media has been treat- 
ed theoretically by Akhiezer,' Kittel and Abra- 
hams,’ and others® in connection with ferromag- 
netic relaxation times. Their work showed that 
in the scattering process of a single phonon with 
a single spin wave (magnon), both energy and line- 
ar momentum were conserved. Turner‘ has ob- 
served that the linewidths of certain spin waves 
in YIG are anomalously large, which he attributes 
to the fact that these spin waves have the same 
frequencies and wavelengths as some phonons, 
and hence interact strongly with these phonons 
and relax more rapidly that other spin waves. 

It is not obvious that the theories which were 
derived for bulk insulating ferromagnets need ap- 
ply to the case of thin films. Near the surface of 
the material there is no translational invariance 
perpendicular to the surface so that linear momen- 
tum need not be conserved in this direction. Thus, 
for films whose thickness is on the order of the 
magnon wavelength or the acoustic wavelength, one 
might expect that it is not necessary that the pho- 
non wavelength equal the magnon wavelength in 
order that interaction occur. To verify these ideas 
experimentally we have shown directly that a se- 
ries of standing spin waves can be excited by pho- 
nons whose frequency equals that of the spin waves, 
but whose wavelength is considerably larger than 
the magnon wavelengths. 

The most direct means of observing spin waves 
is by the method of spin-wave resonance. As ex- 
plained by Kittel® and observed by Seavey and Tan- 
nenwald,°® the pinning of spins at the surfaces of a 
film will allow the excitation of standing spin waves 
by uniform fields. These spin waves precess at a 
frequency given by 


w =y(H-40M)+2Ayk?/M, (1) 


when the external magnetic field H and the wave 
vector k are perpendicular to the film surface. 
Hence y = gyromagnetic ratio, M = saturation mag- 
netization, A = exchange constant, k= magnon wave 
number =p1/L, where L=film thickness and p is 
an odd integer. The experiment (similar in tech- 
nique to one reported by BO6mmel and Dransfeld)’ 
was performed to detect the excitation of these 
spin waves by phonons. 

A film of Permalloy (80% Ni- 20% Fe) of thick- 
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ness 2600 A was evaporated onto the end of an 
AC -cut quartz cylinder. This end of the cylinder 
was placed in the magnetic field region of a rec- 
tangular microwave cavity resonant at 8.9 kMc/sec 
at 4.2°K. The other end of the rod was situated in 
the electric field region of a re-entrant cavity that 
was resonant at the same frequency. First, the 
spin-wave spectrum of the film was measured by 
using the rectangular cavity as part of a micro- 
wave spectrometer. Six strong absorption lines 
were observed. When these lines were fitted by 
the dispersion relation (1) they were found to cor- 
respond to the resonances for p=0, 9, 11, 13, 15, 
and 17. Apparently the spin-wave modes p=1 
through 7 are strongly attenuated, possibly by ed- 
dy current damping.* Then the spectrometer was 
removed and a superheterodyne receiver suitable 
for receiving weak pulses was connected to the 
rectangular cavity. By pulsing the re-entrant cav- 
ity with 10 watts of power of 1- sec pulse width, 
transverse phonons were initiated in the quartz 
bar. These propagated down the bar and entered 
the film which was now in a quiescent cavity. The 
power in the rectangular cavity was observed as 
a function of magnetic field. Power was observed 
at those values of the field at which electromag- 
netically excited spin-wave resonances had previ- 
ously been observed. Thus, the microwave pho- 
nons had excited the magnetization into precession 
in the several spin-wave modes, and at resonance 
the precessing magnetization produced measurable 
excitation of the cavity. 

This demonstrates that phonons of frequency 
8.9 kMc/sec and wavelength ~3400 A interact with 
spin waves of the same frequency and with a num- 
ber of different wavelengths, all less than 1000 A. 
Only four of the spin waves were sufficiently 
strongly excited acoustically to be observed. The 
relative amplitudes of the power output at the 
acoustically excited lines were 1:0.25:0.12:0.05. 
When observed by electromagnetic excitation the 
relative amplitudes of the first four lines were 
1:0.29:0.37:0.13. This indicates that the relative 
phonon excitation is increased as the acoustic and 
magnon wavelengths approach equality. 

The theoretical expectation? is that phonons are 
coupled to spin waves via a magnetoelastic energy 
that is also responsible for magnetostriction. In 
order to investigate this, experiments similar to 
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the one described above were performed on two 
films whose magnetostriction constants, A, were 
considerably different. The average magneto- 
strictive constants of 83Nil7Fe Permalloy® are 
approximately zero, but the average magneto- 
striction of the composition 76Ni24Fe is \ =10 
x10°°, which is greater than that of 83-17 Per- 
malloy by at least a factor of ten. A sample was 
prepared by evaporating a film of 83-17 Permal- 
loy onto one half of the end surface of an AC -cut 
quartz rod. A 76-24 Permalloy film was evapor- 
ated on the other half of the same surface. (This 
was done to assure that the same amount of pho- 
non power entered each film.) The thicknesses of 
the films were made as nearly equal as possible, 
2270 A for the 76-24 Permalloy and 2400 A for 

the 83-17 Permalloy. The electromagnetically 
excited ferromagnetic resonance of this compound 
sample showed two distinct spin-wave spectra cor- 
responding to the two “half-moon” shaped films. 
Because the 47M of the two samples differ by 

more than a thousand gauss the spectra do not 
overlap. The amplitudes and linewidths of the 
main (k=0) lines were approximately the same 

for both compositions when the resonances were 
excited electromagnetically. The phonon-excited 
resonances were quite different, however. It was 
found that spin-wave resonance was excited in the 
76-24 Permalloy (which has the larger magneto- 
striction), but to our experimental sensitivity no 
ferromagnetic resonance was excited in the 83-17 
Permalloy. The phonon spin-wave interaction was 
at least a factor of 15 greater in the 76-24 Per- 
malloy than in the 83-17 Permalloy. It is unlikely 
that the small difference in film thickness has a 
large effect on the excitation; the interaction of 
magnons and phonons seems directly related to 

the magnetostriction. 

From these results some qualitative conclusions 


can be drawn about the relaxation processes in 
Permalloy films. The spin-phonon interaction was 
shown to be strongly dependent on the magneto- 
striction of the films near the 80-20 composition, 
yet the linewidths of the resonances were approxi- 
mately the same (=40 gauss) for all compositions 
studied. This means that spin-phonon processes 
cannot be the major source of the linewidths; 
there are apparently more rapid processes short- 
ening the lifetime of the spins. This lends further 
experimental support to the conclusions of Kittel 
and Abrahams? and Bloembergen and Wang’® that 
direct spin-lattice interactions are not the most 
important source of linewidths in ferromagnetic 
metals at temperatures well below the Curie tem- 
perature. 

It is a pleasure to acknowledge helpful conversa- 
tions with R. Landauer, A. H. Nethercot, Jr., and 
P. Seiden. We are indebted to W. Kately for film 
preparations, and to I. L. Gelles for help with the 
microwave spectroscopy. Particular thanks are 
due J. P. Anderson for his able and enthusiastic 
technical support. 
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SILICON DIVACANCY AND ITS DIRECT PRODUCTION BY ELECTRON IRRADIATION 
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To date two defects produced in radiation dam- 
age of silicon have been identified.’»? These de- 
fects are a vacancy-oxygen pair and a vacancy- 
phosphorous pair. They were identified largely 
by their associated electron spin resonance spec- 
tra and have been labeled the Si-A and Si-E cen- 
ters, respectively. The studies associated with 
these identifications demonstrated that the vacan- 
cy is mobile in a room-temperature irradiation 
experiment and emphasized the importance of 
impurities in influencing the radiation-damage 
processes. In this Letter we report the identifi- 
cation of an intrinsic defect, which is stable at 
room temperature, the divacancy. We also demon- 
strate that this defect can be produced directly in 
a high-energy electron irradiation (that is, the 
migration and agglomeration of vacancies is not 
required for its formation). In addition, a marked 
anisotropy is observed in the production rate with 
respect to the electron beam direction, giving in- 
sight into the detailed production process. 

The divacancy is shown in Fig. 1. It consists of 
two vacant, nearest-neighbor, lattice sites, A 
and B, each surrounded by three occupied lattice 
sites, 1-6. The divacancy identification has been 
made by studies of two associated electron spin 
resonance spectra, previously designated the 
Si-J and the Si-C centers.” These spectra are 
interpreted as arising as follows: The electrons 
on atoms 2 and 3 and on atoms 5 and 6 each pair 
in molecular bonds as shown in the figure. The 
Si-J center is observed in p-type material and is 
interpreted as arising from a single unpaired elec- 
tron in the extended orbital between atoms 1 and 4. 
In this state, the divacancy is a singly ionized do- 
nor. The Si-C center is observed in high-resis- 
tivity silicon and is interpreted as arising from a 
third electron in the antibonding orbital between 
atoms 1 and 4, the two bonding orbitals being filled. 
In this state the divacancy is a singly ionized ac- 
ceptor. The Si-C center is not observed in low- 
resistivity n-type material, indicating that the 
divacancy can also accept an additional electron. 

It thus can act as both a single donor and a double 
acceptor. 

Some of the arguments leading to this identifica- 
tion are summarized briefly as follows: 

1. These spectra have anisotropic g tensors as 
shown in Table I. Both spectra exhibit resolved, 
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FIG, 1. 
labeled A and B represent the two adjacent vacant lat- 


The silicon divacancy. The dashed circles 





tice sites. The electrons associated with the defect 
can be regarded as forming the bonds shown in the 
figure. The axes for the spin-resonance data in Table 
I are indicated in the figure. 


axially symmetric, hyperfine interactions with 
two equivalent silicon atoms. The hyperfine con- 
stants are also shown in Table I. The resonance- 
electron orbital between atoms 1 and 4 (see Fig. 1) 
accounts for the simple (111) hyperfine interac- 
tions with the same axis for two silicon atoms, 
while the g tensor, reflecting the symmetry of the 
wave function as a whole, is tilted away from this 
direction. 

2. For the divacancy shown in Fig. 1, in addi- 


Table I. Spin Hamiltonian constants for the diva- 
cancy spectra (see Fig. 1). 















0 si? hf 
Spectrum g(+0.0003) (See Fig. 1) (107 cm™) 
Si-C 81 =2.0012 
82= 2.0135 6 z}° A, =79 +2 
g3= 2.0150 A,=56 +2 
Si-J  —_g = 2.0004 
£2= 2.0020 7.7+0.5° A, =67.8 +0.5 
g3= 2.0041 A, =40.0+0.5 
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FIG. 2. R vs the average energy of the incident 
electron beam, E. R is the ratio of the number of 
divacancies with the vacancy-vacancy direction along 
the incident beam (111) direction (say the [111] of Fig. 
1) to the number along one of the other three equivalent 
(111) directions (the [11i]). 


tion to the orbital between atoms 1 and 4, as 
shown, equally likely configurations involve the 
orbital being between atoms 2 and 5 and between 
atoms 3 and 6. For both the Si-C and Si-J cen- 
ters, motion between these configurations has 

been detected both by stress-induced alignment 

and by linewidth studies, confirming this predic- 
tion of the model. The activation energy for this 
motion is ~0.06 ev for both centers. 

3. Preferential alignment of the vacancy-vacan- 
cy direction has been quenched in by stressing at 
elevated temperatures. The magnitude of the 
alignment as well as the kinetics of reorientation 
is identical as monitored by both centers. The 
activation energy for this motion is ~1.3 ev. 

4. Both spectra are observed after an in situ 
irradiation at 20°K, showing that they are due to 
a primary defect, which does not require anneal- 
ing for formation. 

That the divacancy is produced as a primary 
process has also been demonstrated in an inform- 
ative way as follows: 

Oriented samples 0.010 in. x0.1 in. x0.6 in. were 
cut from p-type (~3 x10** B/cc) pulled silicon with 
the thin dimension a (110) direction. For conven- 
ience we will consider this thin dimension to be 
the [011] of Fig. 1. These samples were irradiated 
at room temperature with a high-energy electron 
beam incident along the [111] direction, and the 


Si-J center resonance observed. By observing 
the relative intensities of the multiplets of the spin 
resonance spectrum it is possible to count the di- 
vacancies in each of the possible orientations. 
The ratio, R, of the number with the vacancy- 
vacancy axis along the beam direction ({111}) to 
the number in one of the other three equivalent 
(111) directions (i.e., [111]) is shown in Fig. 2 
as a function of bombarding energy. As can be 
seen, there is a substantial anisotropy in the 
production rate, which increases with decreasing 
bombarding energy. 

This behavior can be qualitatively understood in 
terms of a conventional model for divacancy pro- 
duction: A high-energy electron collides with the 
nucleus of an atom and imparts a recoil energy to 
the atom. The recoiling atom collides with a near- 
est neighbor atom, causing it to be displaced into 
the lattice as an interstitial. In addition the initial 
recoil atom must retain enough energy so that it 
too becomes an interstitial. In this way, two adja- 
cent vacancies, a divacancy, are created. The en- 
ergy the primary recoil atom receives in the colli- 
sion with a high-energy electron, and the direction 
of recoil are correlated, being determined by the 
Mott-Rutherford scattering formula. The highest 
energy recoil is along the beam direction. Con- 
sequently as the beam energy is reduced, a point 
is ultimately reached where sufficient recoil en- 
ergy is available to produce a divacancy along the 
beam direction but not in the (111) direction 70° 
32’ away from the beam direction, which suggests 
the tendency seen in Fig. 2. 

There are several considerations which make 
the data in Fig. 2 provisional (e.g., local heating 
of the lattice during the defect production process, 
finite sample thickness, and the fact that the ir- 
radiation was at a temperature where single va- 
cancies could aggregate). However, all of these 
would reduce the observed anisotropy, so that it 
is at least as large as shown in Fig. 2. 

In summary we observe the divacancy as the re- 
sult of the direct production process. It is pro- 
duced at a rate ~5% of the vacancy production rate, 
as monitored by the Si-A center, for a 1.5-Mev 
irradiation,® and plays an important part in room- 
temperature radiation-damage processes. 
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barding energy for the Si-J center indicates a consid- observation, Bemski independently concluded that the 
erably higher threshold energy for divacancy production Si-C center was to be associated with a multiple defect, Lin 
than that for the single vacancy, as is to be expected. suggesting the divacancy as a possibility. — 
I; 
I; 
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BOUND EXCITON COMPLEXES = 
D. G. Thomas and J. J. Hopfield 
Bell Telephone Laboratories, Murray Hill, New Jersey valu 
(Received September 11, 1961) 
Exciton complexes in which excitons are bound This four-level system gives rise to two pairs of 
to neutral or charged donors or acceptors have lines, one with a separation controlled by |g, - gy wher 
been described by Haynes’ and Lampert.? In crys- and the other controlled by |g, +g ,!, where u and H, a 
tals with a direct band gap, such states are ob- l refer to the upper and lower states, respectively. com| 
servable as sharp absorption and emission lines In hexagonal CdS free electrons occupy an S-like stru' 
occurring at wavelengths longer than that of the J=2T, band, and free holes (from the top valence act s 
lowest energy intrinsic exciton state. This Letter band) a P-like J=3 TI, band.* The free electron For 
shows that the Zeeman effect of these impurity has an almost isotropic g value of -1.7,. Because 4 tre 
complex states can be used to identify which type the conduction band is simple, this g value will doub 
of center a given absorption line represents. In be almost independent of the binding of the elec- bitr 
this way, the g values of donors and acceptors, tron in a center. Symmetry requires that the g occu 
the binding energies of the various complexes, ble. 
and even information about the energy levels and ue 0 
values of excited states of donors and acceptors a ' it is 
roe be obtained. The electron-hole spin-spin (ex- @=+|O4- |O-+/0+- Ou 
change) interaction can be measured in some of , a Cds 
these states. The experimental results discussed z S cern 
are for CdS, but the method of investigation seems FE Q Tabl 
to be promising for a wide class of semiconducting 5 & pola 
crystals. Previous work has established the ener- 9% 3 
gies of the intrinsic exciton lines* in CdS, but the ‘\ @- © + @ © 
explanation of the lines occurring at lower ener- NEUTRAL | NEUTRAL | IONIZED | IONIZED 
‘ : ? DONOR | ACCEPTOR | DONOR | ACCEPTOR 
gies has remained unsatisfactory.**® 
Excitons can be bound to neutral donors or ac- HO 
ceptors to form molecular states, or to ionized 5 $9484 vennens gwen 
donors or acceptors to form molecular ion states, (ae 
as illustrated in Fig. 1. The neutral donors and bonds [Su- S| 
acceptors are doublet states (S=4), and this is 
also true for the complexes formed from these 2 is } l9u+9.| Fl 
states since the “bonding” electrons or holes must, ts Heh es and 
have antiparallel spins. (It is well known that for ton 
fluo 
molecular hydrogen, if the “bonding” particles FIG, 1. Schematic diagram illustrating the ground Qe 
have parallel spins an unstable system results. and excited states of the various transitions discussed The 
Precisely the same considerations apply to the in the text. The g values of the states arising from the be : 
exciton complexes discussed here.) The states neutral centers are determined by the g value of the line 


odd particle, and the splittings of these states are as 
illustrated. The ionized centers show no splitting in 
the ground states; the splitting of the excited states of 
these centers with c 1 H is complicated by the presence 
split, as shown in Fig. 1, with a g value equal to of a spin-spin interaction energy. The charges within 4 
that of the odd particle in the particular complex. circle indicate a chemical center. 


considered here are ground states and so have no 
orbital angular momentum. In a magnetic field 
the states derived from the neutral centers will 
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Table I. CdS—1.6°K, 





Energy below 
Energy exciton A(2.5537 ev) 
(ev) (ev) 





0.0177 
0.0066 
0.0038 


2.53595 (4888.5 A) 
2.5471 (4867.15 A) 
2.5499 (4861.7 A) 








value of the hole has the form 
8), =8) i cos6, 


where 6 is the angle between the magnetic field, 
H, and the c axis. The hole, belonging to a rather 
complicated and almost degenerate valence band 
structure, can have a g value sensitive to its ex- 
act state of binding. For the free hole, gp,y = -41.15. 
For cllH selection rules show that the 1-3 and 2- 
4 transitions will occur optically, resulting in a 
doublet split with a g value of |g, -g,|. At an ar- 
bitrary value of 6, however, mixing of the states 
occurs, and a full quartet of lines becomes visi- 
ble. The additional pair is separated with a g val- 
ue of IZ, + &7\ . From measurements of this type 
it is easy to derive the individual g values. 

Our study of “pure” vapor-grown crystals of 
CdS of thickness 10-100 yw has been chiefly con- 
cerned with three prominent sharp lines listed in 
Table I. These lines are strongly active for light 
polarized perpendicular to the c axis, which is al- 


so the polarization of the lowest energy exciton, 
exciton A. At 1.6°K in selected crystals the lines 
may have a width of 0.6 cm™ or less. They can 
be seen at 20°K but not at 77°K. J, and J, are read- 
ily seen in both absorption and fluorescence, while 
I, is observed in absorption. J, is absent in some 
crystals. Sometimes other lines appear near /,, 
the whole group behaving similarly in a magnetic 
field. Additional lines also often occur near J,; 
some of these behave like J, in a magnetic field 
while others behave like J,. Many other lines oc- 
cur, prominent among which are comparatively 
broad lines (about 6 cm™ wide) seen in fluores- 
cence at photon energies below J,; some of these 
correspond to emission taking place with phonon 
creation. 

Lines /, and J, exhibit linear Zeeman splittings 
which can conveniently be studied in fluorescence. 
The splittings of J, and J, as a function of cos@ are 
plotted in Fig. 2 for H =31000 gauss. Quartets 
are seen when 6 has values other than 0° and 90°. 
Within experimental error there is a linear de- 
pendence between the splittings and cos@, so the 
relations 


8, * 8p -8q* 8,4 O88 


must hold, where g,9 is an isotropic electron g 
value equal to -1.7, for J,, and -1.7, for J,. For 
I,, yy equals -1.7; and for J,, g,y equals -2.7,. 
It happens that gp \ = Ze0 for J,. This near equal- 
ity is the reason that for cllH, with the splitting 
determined by g29 - gp} , no Zeeman effect is ob- 
served for this line.® 
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FIG, 2. Splitting of lines J; 
and /, at 31000 gauss as a func- 
tion of cos@ at 1.6°K, as seen in 
fluorescence. Notice that when 
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9*0° or 6 ~90°, quartets occur. 
The zero of energy is taken to 
be at the center of each group of 
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When /, and J, are observed in absorption for 
cllH and cH, the two lines in a pair sometimes 
have different strengths. The differences in 
strength increase as the magnetic splittings in- 
crease or as the temperature decreases, indicat- 
ing that, for the experimental conditions used, 
thermal equilibrium of the spins is being ap- 
proached in the lower states of the complexes. 
The results are as follows: 


I, I, 





No splitting 
observable 


cllH High-energy 
line strong 


Both lines 
equal in strength 


clH High-energy 


line strong 


It is clear that for c1H there is thermalization 
in the lower state of line J,, but that this does not 
occur for line J,. Since ggg=-1.7¢ and gp), =0, J, 
must arise from a neutral donor, and /, from a 
neutral acceptor. The thermalization results for 
cllH are consistent with these conclusions when 
it is remembered that for J,, gy =8eo, and that 
for 1,, |g, !>\geq!l. It is interesting to note that 
the free hole has g=-1.1,, but that in the exciton 
complex @ =+ the value has increased to -1.7, 
and in the neutral acceptor it is -2.7,. 

It follows that the other lines described above, 
which behave magnetically like J, also arise from 
donor states which are, however, chemically 
different from the center responsible for line /,. 
In a similar way, the other lines which behave 
like J, must, as described below, be molecular 
ion states derived from chemically different 
centers. It is an interesting fact that we have 
seen no other lines which behave like J,. This 
may be because there is only one acceptor state 
present in appreciable quantity, but it seems 
more probable that this is the only “shallow” 
acceptor state present; deeper lying acceptor 
states probably give rise to much broader lines 
as a result of stronger coupling to the acoustic 
phonons of the lattice. 

The fluorescent intensities of J, and J, show 
that spin thermalization is not occurring in the 
upper states, presumably because of the short 
lifetime of these states. 

CdS is normally m type, so acceptors are 
therefore usually compensated (ionized). On 
illumination with band gap light, however, free 
holes and electrons are made. These can be 
trapped at low temperatures by the ionized ac- 
ceptor and donors, respectively, to form neu- 
tral centers on which the exciton complexes 
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can be formed. It is found that when the lines 


are observed in absorption using as little band 
gap light as possible, infrared illumination of 
the crystal causes /, and /, to diminish in strength 
(and in some cases to vanish), while J, increases 
markedly in strength. Since the infrared illumi- 
nation is of such a wavelength as to ionize the ac- 
ceptor preferentially, we may conclude that this 
occurs and that the free hole then ionizes the 
donors. Thus, the neutral donor and acceptor 
concentrations are decreased while those of the 
ionized states are increased. We may there- 
fore conclude that J, corresponds to the forma- 
tion of a molecular ion state. 

The magnetic behavior of J,, although not sim- 
ple, appears consistent with its assignment as a 
molecular ion state. For example, no thermal- 
ization is observed for absorption, the expected 
situation for a singlet ground state. In a magnetic 
field with c1H, J, splits unsymmetrically in the 
sense that the lines do not shift linearly with H; 
at low fields only the high-energy member of the 
doublet is visible, but as the field increases the 
low-energy line is seen with steadily increasing 
strength. Extrapolation of the positions of the 
lines to H =0 shows that there is a zero-field 
splitting of 3.1x10~* ev. Symmetry arguments 
show that with c1H, states with antiparallel 
spins (1S [, states) are mixed with states with 
parallel spins (1S [, states). The zero-field 
splitting is therefore ascribed to a spin-spin 
interaction energy, and since the transition to 
the 1S [, state is forbidden, it is clear that this 
state has an energy lower than the 1S I, state. 
Such an interaction can occur in absorption lines 
due to excitons bound to ionized donors or accep- 
tors, where the complex involves an odd number 
of electrons and of holes. It cannot occur in ab- 
sorption lines due to excitons bound to neutral 
donors or acceptors, since the spins of the two 
like particles will be antiparallel. The spin- 
spin interaction energy of 3.1x10~* ev is sur- 
prisingly large. We have, however, evidence 
from the effect of strain on the intrinsic exciton 
lines that the free exciton has a similarly large 
value. 

The methods described here could be of help 
in determining the types of impurities present 
in many 2-6 semiconductors. Information is also 
derived about the donors and acceptors, as il- 
lustrated here by the determination of the g value 
of the hole on the acceptor (a value that would 
be hard to obtain by paramagnetic resonance 
since Am =3 for a spin transition in the J= 
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band). The dethermalization of the lower state 
population by microwave pumping should be 

easily detected optically. Donor and perhaps 
acceptor g values can therefore be measured 

with microwave precision on crystals too small 

to investigate with normal microwave techniques. 
Finally, there is evidence in the fluorescent 
spectra for transitions which leave neutral donors 
or acceptors in excited states. Such transitions 
provide a means of measuring donor and accep- 


tor excited state energies and magnetic split- 


tings. 





1J. R. Haynes, Phys. Rev. Letters 4, 361 (1960). 

2M. A. Lampert, Phys. Rev. Letters 1, 450 (1958). 

3J. J. Hopfield and D. G. Thomas, Phys. Rev. 122, 
35 (1961). 

4{L. S. Pedrotti and D. C. Reynolds, Phys. Rev. 120, 
1664 (1960). 

5R. G. Wheeler and J. O. Dimmock, Phys. Rev. 
Letters 3, 372 (1959). 





M-SHELL, RELATIVISTIC ELECTRON WAVE FUNCTIONS* 


H. R. Brewer, D. S. Harmer, and D. H. Hay 


Georgia Institute of Technology, Atlanta, Georgia 
(Received August 3, 1961) 


The effects of screening and finite nuclear size 
on M-shell electron wave functions have been cal- 
culated numerically. The ratio of the corrected 
M-shell wave functions to the corresponding point- 
nuclear values (evaluated at the finite nuclear ra- 
dius) is presented here in Fig. 1 for their applica- 
tion in evaluating M -capture transition probabili- 
ties.’ In the past, these orbital capture transi- 
tions were corrected empirically for screening. 
The results presented here indicate the necessity 
for including such corrections and, in addition, 
indicate wide variations between the correction 
factors for different M subshells. 

The correction factors were derived from wave 
functions which were calculated as part of a pro- 
gram for obtaining M -shell internal conversion 
coefficients. The wave functions are available in 
tabular form’ for points outside the nuclear sur- 
face for atomic numbers 55 through 90 in steps 
of 5. 

For the evaluation of these wave functions, a 
uniform charge distribution was assumed inside 
the nucleus, and power series solutions* were 
obtained for f and g (small and large component) 
of Dirac’s equations. These series solutions were 
evaluated at the nuclear radius® to give starting 
values for a numerical solution of the equations 
in the external region. Here the potential due to 
the nucleus and orbital electrons was represented 
by a Thomas-Fermi-Dirac (TFD) potential® and is 
given by 


V=aZg/r+D. 


In this equation @ is the fine structure constant, 
? is the screening function taken from Umeda’s 


tables,” D is the Dirac exchange term,® Z is the 





atomic number, and ¢ is the radial distance. The 
outer boundary condition for f(a) and g(a) is given 
by 


f(a)/g(a) = -[(1 - W)/(1 + W)), 


where W is the trial energy eigenvalue and a is 
the value of y for which g vanishes. Trial values 
of W were varied until the solutions matched both 
the boundary condition at y=a and that at the nu- 
clear surface, and possessed the correct number 
of nodes in the interval. 
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FIG. 1. f, /f, (c) and g, /g, (c) are the ratios of the 
screened, finite nuclear size functions to the corre- 
sponding analytic point-nuclear functions’ evaluated at 
the nuclear radius and plotted as a function of the nu- 
clear charge Z. The numbers tc the right of each curve 
give the value of « as defined in the text. 
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Table I, Energy levels for M-shell electrons computed with finite size and TFD screening corrections compared 


with experimental and Coulomb point-nuclear values. 











(Subshell) W (mc?) 


(1-W) (1-W) (1-W) 
Calculated Experimental® Coulomb point 
(ev) (ev) nucleus (ev) 





0.998542 
0.998511 
0.998014 
0.997868 
0.997612 


0.992863 
0.992501 
0.991493 
0.989854 
0.989251 


745 720 4594 
760 738 

1014 

1089 

1220 


3646 
3831 
4346 








a 
See reference 9. 


Figure 1 shows the ratio of screened, finite- 
nuclear-size f and g to their corresponding Cou- 
lomb, point-nuclear values® at the nuclear radius. 
The calculation errors (truncation and round-off) 
in f and g are believed to be less than 05%. These 
were estimated by inserting a Coulomb potential 
in place of the screened one and beginning the nu- 
merical integration at the nucleus with the Cou- 
lomb, point-nuclear starting values. The results 
were then compared with the analytic solutions® 
to obtain the error estimate above. 

A comparison of the TFD energy levels with the 
Coulomb and experimental values® is given in Ta- 
ble I for the highest and lowest values of Z. The 
deviations between computed and experimental 
values for intermediate Z lay in between those 
given in Table I. The symbol x used to designate 
the levels has the conventional meaning!: 


x = -(1+1) when j=1+3, 


k= 1 when j=1-3. 


The correspondence of the experimental and com- 
puted energy levels over this range of Z lends 
confidence to the use of the TFD screened poten- 
tial. 

In conclusion the authors wish to thank Dr. 
M. E. Rose for suggesting the problem and for 


his kind assistance during the course of the work. 
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ELECTRON-ION RECOMBINATION COEFFICIENTS IN NITROGEN AND IN OXY GEN" 


W. H. Kasner, W. A. Rogers,! and M. A. Biondit 


Westinghouse Research Laboratories, Pittsburgh, Pennsylvania 
(Received September 22, 1961) 


Currently, there is considerable interest in 
predicting free-electron lifetimes in various re- 
gions of the ionosphere following ionizing events, 
such as solar flares. An important process in 
determining the electron loss is recombination 
between electrons and the positive ions present 
in the upper atmosphere. It has been reasonably 
well established from analysis of ionospheric ob- 
servations? and from laboratory experiments? ® 
that the recombination coefficient a (defined by 
dn, /dt= -an,n,, where no and n, are the electron 
and positive ion densities, respectively) may be 
as large as 10~® to 10~* cm*/sec and that the cap- 
ture process is dissociative recombination.‘ 

Since no adequate theoretical calculations of the 
rate of the dissociative recombination reaction, ° 
e.g., 


(XY) +e——-» (XY) " x"+¥, (1) 


—_—_— 
unstable 
are available, appeal is generally made to lab- 
oratory measurements of the recombination co- 
efficients in atmospheric gases.*»* Unfortunately, 
previous measurements were made without identi- 
fication of the ions responsible for the recombina- 
tion and under conditions where more complex ion 
types, e.g., O,* in oxygen and N, and N,* in ni- 
trogen, are observed. This Letter presents the 
preliminary results of an experiment involving 
microwave determinations of the decay of elec- 
tron density during the afterglow following ion- 
ization of gas in a cavity,” coupled with the simul- 
taneous identification, by an rf-type mass spec- 
trometer, of the afterglow positive ion currents 
to a wall of the cavity. From these measurements, 
recombination coefficients appropriate to particu- 
lar positive ion species, e.g., N,*, can be as- 
signed, 

The apparatus consists of a rectangular, stain- 
less steel, microwave cavity containing in one 
wall a quartz iris for coupling microwave energy, 
in the opposite wall a 0.008-inch diameter effusion 
orifice leading to a Boyd-type® rf mass spectrome- 
ter and its differential pumping system, in a third 
wall an optical viewing port, and in a fourth wall 
a lead to an ultra-high vacuum, gas-handling sys- 
tem. The gas-handling system, the cavity, and 
the mass spectrometer differential pumping sys- 
tem are baked at 350°C for 14 hours to minimize 
impurity effects. The high-purity gas samples, 
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FIG, 1. Electron and N.” ion decay in the afterglow 
of N»-He microwave discharges having constant nitrogen 
pressure, 0.006 mm Hg, and variable helium pressure. 
The N,* wall currents are normalized to the respective 
electron density decay curves at an afterglow time of 
approximately 2 milliseconds. 


contained in steel pressure tanks,° are fed into 
the system through controlled leaks. 

In order to carry out measurements at low con- 
centrations of nitrogen or oxygen, an inert buffer 
gas is used to inhibit particle loss by diffusion to 
the walls. It was observed that in either pure 
nitrogen or nitrogen-helium mixtures, the pre- 
dominant afterglow ions for nitrogen pressures 
in the range 0.1 to 7.0 mm Hg were N,* and N,*, 
while at nitrogen pressures less than 0.01 mm Hg, 
N,* was the only significant ion in the afterglow.’° 
Similarly, in oxygen, at partial pressures of less 
than 0.005 mm Hg, O,* predominated. 

At low total pressures, see Fig. 1, the pre- 
dominant loss of electrons and N,* ions in the late 
afterglow is by ambipolar diffusion to the walls, 
as indicated by the fact that the rate of particle 
loss varies inversely with gas pressure and ex- 
ponentially with time. The curvature in the top 
curve shows significant electron-ion recombina- 
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tion in the early afterglow. From the known fun- 
damental diffusion length of the cavity, A=1.126 
cm, a Dap value of approximately 900 (cm/sec) 
x (mm Hg) is obtained for N,* ions and electrons 
moving in helium. At this low partial pressure 
of nitrogen, 0.006 mm Hg, N,* and N,* ions are 
essentially absent. The decay of the normalized 
wall current of N,* ions (which, for fundamental 
mode diffusion, is proportional to the average 
density of N,* inside the cavity) follows the aver- 
age electron density, within the accuracy of the 
present measurements." 

In order to study recombination under conditions 
where diffusion loss is unimportant, it is neces- 
sary to increase the pressure of the buffer gas to 
reduce the diffusion rate. Neon is more suitable 
than helium, the gas conventionally used, ®»!” since, 
at a given pressure, neon inhibits diffusion more 
strongly"* and effuses into the mass spectrometer 
more slowly than helium, thus reducing choking 
of the mass spectrometer pumping system. The 
ionization potential of neon is sufficiently high 
that one observes only ions of the minority nitro- 
gen or oxygen constituent. 

Measurements of electron-ion recombination 
have been carried out under conditions where one 
particular ion species (except for the N,*-N,* 
combination) predominates during the afterglow. 
Under these conditions, and when diffusion loss 
can be neglected, the solution of the electron den- 
sity equation becomes (1/n,)=(1/ng9)+ at. Ex- 
amples of data obtained for recombination loss of 
N,* ions and electrons and for O,* ions and elec- 
trons are shown in Fig. 2. The neon gas pressure 
is sufficiently high so that the ratio of recombin- 
ation loss to fundamental mode diffusion loss is 
large over the entire range, varying from ~200 
to ~10 in the data shown. From the initial straight 
portions of the curves, covering a factor of ~10 
in electron density, taken at several neon partial 
pressures, we obtain the values a(N,*)=(5.9+ 1) 
x107? cm’/sec and a(O,*)=(3.8+1)x10~" cm/sec. 
The uncertainty in these values arises largely 
from imperfect knowledge of (a) the correction of 
the recombination curve for diffusion effects in 
this geometry, * (b) the spatial density distribu- 
tion of the electrons, which affects the determina- 
tion of absolute electron density from the meas- 
ured cavity frequency shifts,’ and (c) the time-re- 
solved afterglow ion currents, which indicate the 
importance of any other positive ions in the after- 
glow. Electron attachment in oxygen should not 
be a significant factor in the recombination analy- 
sis, since the oxygen partial pressure was so low, 
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FIG, 2. Time variation of the reciprocal electron 
density in the afterglow of high-pressure N,-Ne and 
O,-Ne microwave discharges. 


i.e., less than 107° mm Hg. 

At higher partial pressures of nitrogen, N,* and 
N,* ions appear in nearly equal abundance and 
seem to vary together with both gas pressure and 
time during the afterglow. The curves of 1 /Ne vs 
time are linear in this case also, indicating an ef- 
fective recombination coefficient of the order of 
2x107* cm’/sec. 

The present results indicate that at the pressures 
used in the earlier work, ?»® the ions under study 
were not N,* and O,*, as had been assumed. Our 
preliminary value for a(N,*-N,*) agrees with that 
reported in the earlier nitrogen work? and with 
some recent results of Van Lint et al. ,’* while our 
value for a(N,") is substantially smaller, thus in- 
dicating the need for ion identification in atmos- 
pheric gas studies of this type. Studies of the re- 
combination coefficients of other ions, such as 
(NO)*, of importance in upper atmosphere prob- 
lems, are being considered. 
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tPresent address: Physics Department, Thiel Col- 
ap Greenville, Pennsylvania. 

University of Pittsburgh and Westinghouse Research 
Laboratories, Pittsburgh, Pennsylvania. 

‘See, for example, D. R. Bates and H. S. W. Massey; 
Proc. Roy. Soc. (London) A187, 261 (1946); H. S. W. 
Massey, Advances in Physics, edited by N. F. Mott 
(Taylor and Francis, Ltd., London, 1952), Vol. 1, 
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PHOTOPRODUCTION OF SINGLE NEUTRAL PIONS FROM HYDROGEN AT 60° 
IN THE 600-1100 Mev REGION 


R. Diebold, R. Gomez, R. Talman, and R. L. Walker 
California Institute of Technology, Pasadena, California 
(Received September 20, 1961) 


In a recent experiment Cortellessa and Reale’ 
measured the differential cross section for the 
process y+p-~7°+p in the neighborhood of the 
second resonance at a center-of-mass angle for 
the neutral pion of 57°. Their results, obtained 
with a proton telescope in coincidence with a lead 
glass Cerenkov counter used to detect one of the 
7 decay photons, are considerably different from 
those which had been reported previously. Earlier 
data were obtained by Vette” using a proton mag- 
netic spectrometer, by Stein and Rogers® and Wor- 
lock* using proton telescopes, and by DeWire, 
Jackson, and Littauer® using a proton telescope 
plus Cerenkov gamma-ray counter. Some of the 
difference can be attributed to the better resolu- 
tion of Cortellessa and Reale which was +30 Mev 
in incident photon energy, about half that of the 
previous experiments. However, even taking in- 
to account this difference in resolution, a serious 
disagreement exists between the experiment of 
Cortellessa and Reale and the earlier ones in that 
Cortellessa and Reale obtain values for the cross 
section which are, in general, considerably lower 





than those previously obtained. Furthermore, 
Cortellessa and Reale find a sharp peak in the 
cross section at a photon energy 700 Mev, which 
is about 50 Mev lower than the energy correspond- 
ing to the “second resonance” peak observed in 
m-nucleon scattering.® Although the excitation 
curve at any given angle need not have a maximum 
at the same energy as the total cross section, the 
angular distributions for 7° photoproduction do not 
seem to be changing rapidly in this region, so that 
one might expect the peak at 57° to be character- 
istic of the total cross section and to be centered 
closer to 750 Mev than 700 Mev as measured by 
Cortellessa and Reale. (The fact that 7* photo- 
production shows a peak at 700 Mev is probably 
the result of an interference effect with the meson 
current term which does not appear in 7° photo- 
production.”) 

Intrigued by this apparent discrepancy and by 
the interesting shape of the curve reported by 
Cortellessa and Reale, we have measured the 7° 
cross section in this region and extended the meas- 
urements through the third resonance using the 
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bremsstrahlung beam of the Caltech Synchrotron. 
The measurements were made by counting recoil 
protons alone and also by counting proton—gamma- 
ray coincidences. The protons were detected by a 
wedge-shaped, uniform field, magnetic spectrom- 
eter and associated counters.® They were identi- 
fied both by a time-of-flight measurement and by 
the pulse heights in three scintillation counters 
located at the magnet focus. Decay photons from 
the 7s were detected in two different ways. One 
was by a large lead glass Cerenkov counter (30 
cm high, 36 cm wide, and 36 cm deep) placed in 
the 7° flight direction with the front edge 65 cm 
from the liquid hydrogen target. The second was 
by using a 1.18-cm lead converter and a scintilla- 
tion counter (16 cm x20 cm) placed 60 cm from 
the target directly in front of the Cerenkov counter. 
Coincidences between the scintillator and the Cer- 
enkov counter were then used to measure gamma 
rays in the smaller solid angle subtended by the 
lead converter. An anticoincidence counter was 
not used to eliminate charged particles since it 
was found that such a counter accidentally vetoed 
some 7° events because of the numerous low-ener- 
gy charged particles present. A separate meas- 
urement was made of the number of “gamma-ray” 
counts produced by charged particles. We did not 
limit the active part of the Cerenkov counter with 
an absorbing aperture since such an aperture 
would not have been well defined without an anti- 
coincidence counter. Since the counter subtended 
such a large solid angle, very few of the 7°-decay 
y rays passed near its edge, and edge effects were 
thus quite small. 

For the backward proton center-of-mass angle 
under investigation, it was difficult to obtain good 
photon energy resolution. Although we reduced 
the proton momentum acceptance interval to half 
that used by Vette, the resolution was decreased 
by only 20% because the angular acceptance of 
the magnet (+1°) contributed a great deal to the 
spreading of the resolution. Sizable contributions 
were also made by multiple scattering in the tar- 
get and by slowing down in the hydrogen. Typical 
values of the resolution, R(®) (full width at half 
maximum), of the experiment are: R(585)=70, 
R(780) = 83, R(1073)=112 Mev. 

Three different counting rates were monitored 
and each was used separately to calculate the 
cross section: (1) the recoil protons alone; (2) 
protons in coincidence with the Cerenkov counter; 
(3) protons in coincidence with both the Cerenkov 
counter and the scintillation counter in front of 
the Cerenkov counter. The geometric efficiencies 
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of the y-ray counters were calculated using a Mop. 
te Carlo program. The Cerenkov counter was Cal- 
culated to be 90% efficient at k= 585 Mev, in- 
creasing to 100% atk=731. The geometric effj- 
ciency of the scintillation counter ranged from 
52% at k=585 Mev to 92% at R=926. The con- 
version efficiency of the 1.18 cm of lead in front 
of the scintillation counter was estimated to be 
80 %. 

Empty-target backgrounds of protons alone 
were approximately 30% of the full target count- 
ing rate. However, the Cerenkov counter coin- 
cidence requirement reduced this background to 
about 4%, and thus improved the counting statis- 
tics by a factor of about 1.5. Requiring both the 
scintillation and Cerenkov counters led to poorer 
statistics due to the lower efficiency for detecting 
the 7° decay photons. 

The set of points shown in Fig. 1 was run by 
varying the bremsstrahlung end point, E,, with 
fixed kinematic settings of the spectrometer and 
counters corresponding to a photon energy, k,- 
of 800 Mev for single 7° photoproduction. This 
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FIG. 1. Number of counts per equivalent quantum as 
a function of bremsstrahlung end point, Ey, with fixed 
kinematic settings of the spectrometer and counters 
corresponding to k= 800 Mev and @=60° c.m. for the 
reaction y+p—7°+p. Folding the resolution given by 
the Monte Carlo program together with the brems- 
strahlung spectrum as measured by J. Boyden [Ph.D. 
thesis, California Institute of Technology, 1961 (un- 
published)] and the cross section as measured by this 
experiment gives the solid line. The errors shown are 
those of the counting statistics only. For a given Eo 
the errors for the three methods used are not independ- 
ent. 
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gave a check on the resolution, on the correlation 
petween k as determined by the magnet system 

and E, as determined by magnetic field measure- 
ments of the synchrotron, and on the contribution 
from pion pairs. A comparison of the proton plus 
Cerenkov counter data with the calculated curve 

of Fig. 1 shows a discrepancy between k and E, 

of 17+ 3 Mev. Measurements at other points have 
indicated a similar displacement. This small dis- 
crepancy is not yet understood. The rate of rise 
of the experimental points for E, < 900 Mev agrees 
quite closely with that of the curve verifying the 
resolution calculations. The rise above E,=900 

is presumably due to pion pairs. Near the thresh- 
old for pion pairs (about 950 Mev for the equip- 
ment set at k= 800 Mev) the 7’s come out in such 

a direction that most of them pass through the 
Cerenkov counter. To measure the charged parti- 
cle contribution to the counting rate at higher E,, 
the front scintillation counter pulse heights were 
examined at E,=1200 Mev and also at 900 Mev 
where it was assumed that no charged particles 
were counted. Charged particles gave a minimum- 
ionizing pulse whereas the showers produced in 
the lead converter by the 7° decay photons in gen- 
eral gave at least twice-minimum-ionizing pulses. 
Subtracting the charged particle contribution at 
1200 Mev estimated in this way gave the point 
shown in Fig. 1, indicating that within the statisti- 
cal accuracy of this measurement, all of the rise 
above 900 Mev may be ascribed to charged parti- 
cles, presumably pions, counting in the Cerenkov 
counter. 

Figure 1 also shows the extent to which cross 
sections obtained from the three different count- 
ing rates described above agree with each other. 
The counting rates of protons alone, properly 
normalized, were consistently higher than those 
with a gamma-ray coincidence. Furthermore, at 
£,=700 Mev, which is below the limit of the reso- 
lution function, the proton counting rate (with emp- 
ty target background subtracted) did not vanish. 
This indicates a background from hydrogen pro- 
duced by lower energy photons. A mass-spectro- 
graph analysis of the hydrogen showed no impurity 
which could account for the background. The 
mechanism producing these protons is not fully 
understood. Some of the protons probably come 
from particles, photoproduced by the beam, which 
Scatter in the hydrogen (for example, y+p+1*+n 
followed by n+p+n+p). Calculations are being 
made to determine whether such scattering can 
completely explain the effect. The proton plus 








converted j-ray results tend to be lower than the 
proton plus Cerenkov results; more detailed cal- 
culations of the lead conversion efficiency may 
clear up this slight discrepancy. For these rea- 
sons, the proton plus Cerenkov data are consid- 
ered to be the most precise and accurate of the 
three sets of data. 

The results of the experiment are shown in Fig. 
2. Assuming a systematic error in the proton 
magnetic spectrometer settings of the magnitude 
indicated by the discrepancy shown in Fig. 1, the 
points were plotted at energies of approximately 
2% lower than those given by the magnet settings. 
Empty-target backgrounds have been subtracted, 
and the data have been corrected for the absorp- 
tion and scattering of the protons in the magnet 
system and for the Cerenkov geometric efficiency. 
Pulse-height spectra from the Cerenkov counter 
were measured and indicated no loss of efficiency 
from the electronic bias. The experiment was run 
with E,=k+100 Mev, so that pion pairs were kine- 
matically excluded and no correction was neces- 
sary. Due to the similarity of kinematics, elas- 
tic scattering of photons by protons (proton Comp- 
ton effect) could not be excluded. Measurements® 
have shown this effect to be of the order of 3% of 
the single 7° production. For 585<k <878 Mev 
the points in Fig. 2 were obtained from the proton 
plus Cerenkov data. Above 900 Mev the Cerenkov 
counter was not used and the results were obtained 
using protons only. As a correction for the proton 
background from lower energy photons and proc- 
esses other than 7° production, a fixed amount, 
0.31 yb/sr, was subtracted from these results 
(0.31 yb/sr is the average difference between the 
results from protons alone and the proton plus 
Cerenkov results; this difference was remark- 
ably constant over the range in which the Cerenkov 
counter was used). Our results, while agreeing 
substantially with the early experiments?“ as 
typically represented by that of DeWire et al., 
present a much clearer picture of the second 
resonance. These results disagree with those of 
Cortellessa and Reale. The third resonance ap- 
pears to be much less prominent than the second 
resonance and has a width comparable with that 
of the second resonance. The widths are not well 
determined because of our wide energy resolution. 

Using protons alone, data are being obtained at 
c.m. angles 90° and 120°. The results obtained 
thus far show a shape at 90° similar to the one at 
60°, largely agreeing with the values obtained by 
the early experiments once the difference in reso- 
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EVIDENCE CONFIRMING THE T =0 THREE-PION RESONANCE* 


Nguyen Huu Xuong? and Gerald R. Lynch 


Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received September 15, 1961) 


We present here some evidence for a T =0, 
three-pion resonance with data from events of 


the type 
p+p>3nt+3n~ +7°, (1) 


produced by antiprotons of 1.61 Bev/c in the 72- 
in. hydrogen bubble chamber.’ The energy and 
the width of the resonance agree very well with 
the w meson found by Maglié et al. who analyzed 
the 27*+27~+7° annihilations in the same experi- 
ment.” 

We have analyzed 551 six-pronged events which 
are a very pure sample of six-prong pion annihi- 
lations: p+p+37*++3n~+nn°. Of these events, 
<1% are 17> +p+p+32~+2n*+n7°, Likewise, Kalb- 
fleisch has shown that <1% of these events are an- 
nihilations involving K mesons.* When these six- 
pronged events were fitted by using the kinemat- 
ics program KICK, 229 were found to have a x? 
less than 5.1 for reactions of type (1) and did not 
fit the hypothesis that no 7° was produced. (The 
latter hypothesis was rejected when the x? for it 
was greater than 30.) An examination of the mis- 
sing-mass distribution as well as the x” distribu- 
tion convinces us that these 229 events contain a 
contamination of less than 5% of the eight-body 
annihilations and less than 3% of the six-body an- 
nihilations. On the other hand, roughly 10% of 
the true seven-body events have been excluded 
from the sample and are grouped with the six- 
and eight-body events. 

We have evaluated the three-body effective 
mass, 


M,= [(E, +E,+E,)? - @, +b, +,)?}”, 


for each pion triplet in the 229 events. Each of 
these events yielded 33 such quantities corre- 
sponding to the charge states: 


1Q|=0: m*1~n° (229x9 combinations), 
1Q| =1: 2*a*a* (22918 combinations), 
1Q| =2: 


We calculated for each value of M, an uncertainty 
6M, by using the variance-covariance matrix of 
the fitted-track variables, which is evaluated by 
KICK. The half-width, T'yego)./2, of the resolu- 
tion function of M, is 9.0 Mev. However, because 


of systematic errors known to exist in our track 


mt*q+7° (229 x6 combinations). 


Number of pion triplets 


reconstruction, our estimate of T'yego),/2 prob- 
ably should be increased by V2 to 'yego1,/2 =13 
Mev. 

Figure 1 is the ideogram of the M, distribution 
for the 229 events. Distributions 1(a) and 1(b) 
are for the charge combinations |Q| =1 and 2, 
respectively. To show the difference between 
the neutral M, distribution and that for |Q| 21, 
we have replotted at the bottom of Fig. 1 both 
the neutral distribution and 9/24 of the sum of 
the |Q| =1 and |Q| =2 distributions. This latter 
distribution we use as an estimate of the phase- 
space distribution. The neutral distribution, 1(c), 
shows a peak at 780 Mev that contains 76 pion 
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FIG. 1. Ideograms of the distributions of the effec- 
tive masses (Ms) of pion triplets for the 229 p+ p— 32* 
+32~+7° events; (a) is for the distribution for the trip- 
lets with |Q|=1, (b) is for triplets with |Q|=2, and 
(c) is for the triplets with Q=0. In (d) the combined 
distributions of (a) and (b) (shaded area) are compared 
with the (c) distribution. The same smooth curve has 
been drawn on (a), (b), and (c). 
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FIG, 2. The density of the population of events on 
the Dalitz plot for five regions of approximately equal 
area as a function of the average distance of this region 
from the center of the Dalitz plot. Both scales are in 
arbitrary units. (a) Events in the energy region of the 
w meson and in a control region. (b) Difference be- 
tween the two sets of points in (a). The three curves 
correspond to the predictions for this distribution on 
the assumptions that the particle has spin and parity 
1~, 1*, or 0. 


triplets above a phase-space estimate of 246. 
Therefore, we confirm the determination of the 
isotopic spin of the resonance as zero. The peak 
has a half-width of about 18 Mev in agreement 
with Maglié et al.* We must conclude, as they 
did, that the half-width of the experimental peak 
is so close to our resolution that the true width 
of the peak is less than 18 Mev and could be zero. 
Of the 229 seven-pion annihilations, 76+ 18 or 
(33 + 8) % of these proceeded via the reaction p 
+p+w+2n*+2n~, for which the cross section 

is 0.6+0.15 mb. 

To determine the spin and parity of this reso- 
nance (or unstable particle), we use the method 
devised by Stevenson et al.,° who have made such 
a determination and conclude that the w is a vec- 
tor meson. We use events of the center of the 
peak (760 to 800 Mev) and events in a control re- 
gion (822 to 878 Mev) of the same charge state. 
About 27% of the events in the peak region belong 
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to the resonance. In Fig. 2(a), we plot the num- 
ber of events per unit area of the Dalitz plot for 
the peak region and for the control region versus 
the distance from the center of the plot. The num- 
ber of events per unit area of the control region 
is normalized to be equal to the estimated back- 
ground in the peak region. Figure 2(b) shows the 
difference of the two sets of points in Fig. 2(a). 
We assume that this difference represents the 
number of events per unit area of the resonance.® 





These data agree very well with the curve pre- 
dicted by a matrix element of a vector meson (1-) 
and not at all with the prediction of a pseudoscalar 
(0-) or axial vector meson (1+). This result of the 
spin and parity (1-) also agrees with that found by 
Stevenson et al.° 

With the spin and parity (1-), we conclude that 
the T =0 three-pion resonance can be the particle 
predicted by Nambu to explain the electromag- 
netic form factors of the proton and neutron,’ 
and also expected in the vector-meson theory of 
Sakurai® or as a member of an octet of vector 
mesons, according to the unitary symmetry the- 
ory.® Chew has pointed out that such a vector 
meson can exist on dynamical grounds as a bound 
state.?° 

We want to thank Professor Luis Alvarez for 
his constant encouragement during the experi- 
ment and for his foresight and effort that made 
the experiment possible. We want also to thank 
Dr. M. Gell-Mann, Dr. B. C. Magli¢, Dr. A. H. 
Rosenfeld, and Dr. M. L. Stevenson for many 
helpful discussions and advice. We wish to ac- 
knowledge the active participation of Dr. C. Tate 
and J. Requa in the analysis of the data. The film 
used in this measurement was obtained in collab- 
oration with Dr. J. Button, Dr. P. Eberhard, Dr. 
G. R. Kalbfleisch, Dr. J. Lannutti, Dr. B. C. Mag- 
li¢é, Dr. M. Pripstein, and Dr. M. L. Stevenson; 
this experiment would not have been possible with- 
out their help. 





*Work done under the auspices of the U. S. Atomic 
Energy Commission. 

Ton leave from the Viet-Nam Atomic Energy Office, 
Saigon, South Viet-Nam. 
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‘The resolution for an ideogram is ¥2 larger than the 
resolution for a histogram. Thus the resolution that 
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